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In this paper, we extend Simon’s criterion for Gaussian states to the multi-mode Gaussian states 
using the Marchenko-Pastur theorem. We show that the Marchenko-Pastur theorem from random 
matrix theory as necessary and sufficient condition for separability. 


The most important question in the quantum informa¬ 
tion is how does one say that the given state is entangled 
or separable? This question has been be answered by 
peres [l| in the finite dimensional case. Horodecki [H has 
shown that this condition is necessary and sufficient for 
separability in the 2x2 and 2x3 dimensional cases, 
but fails for the higher dimensions. In the literature, this 
condition is known as Peres-Horodecki criterion. The 
Peres-Horodecki criterion for the covariance matrix us¬ 
ing Wigner distribution has been given for the Gaussian 
states by Simon Q . It has been shown in the same paper 
that the transpose operation T, which takes every p to 
its transpose fF, is equivalent to a mirror reflection in 
phase space: 

P —> F W{q,p) —)> W{q,-p). (1) 

In this paper we extend the Simon’s criterion to the 
ensemble of n-bipartite system of two modes using the 
Merchnko Pasture law in random matrix theory. The 
Marchenko-Pastur law appears in the famous Marchenko- 
Pastur theorem [^. 

In random matrix theory, the dynamics of the ensemble 
of an infinite dimensional random matrix is described by 
the probability distribution function Q 

V{\i,...,\N)dK = Cne-^^ dk ( 2 ) 

where the matrix H is brought into the following form 
H = E(Aj) - - ^j\ here V{\i) is the 

potential and the Xi are the eigenvalues with i being a 
free index running from 1, 2..A^, then |Ai — Xj\ is the Van¬ 
dermonde determinant, dk = dXi... dXjsi, Cn is constant 
of proportionality and the index /3 = 1,2,4 characterises 
the real parameters of the symmetry class of orthogo¬ 
nal, unitary and symplectic respectively for the random 
matrix H. 
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The invariant ensembles in random matrix theory is 
classified into three class, the Gaussian ensembles these 
matrices are known as the Wigner matrices, the Wishart 
matrices and the two Wishart matrices. For details ref 
0. The Marcenko Pastur’s Quarter-Circle Law is given 
for the Wishart matrices and they are defined as follows. 

Namely, let Xi, • • •, X„ be independent and identically 
distributed random column vectors of 7?.™, with mean 0 
and covariance matrix Im- Let us consider the m x m 
empirical covariance matrix 

1 " 1 

S = - y XiXj’ = where F = (Vi • • • V„). (3) 

n ^' n 


Then the distribution of the eigenvalues of E are con¬ 
sidered under the different assumptions on the number 
of variables m and the number of observations n or m 
is the sample size and n is the dimension of the vectors. 
Their joint probability distribution function of eigenval¬ 
ues is given by the Laguerre polynomials 0. We have 
E{XyY) = Im and the strong law of large numbers in¬ 
dicates that with probability one, 

linin^oo yy = Im- (4) 

n 

In the case m grows with n, here m is the number 
of random variables and the n is the dimension of the 
Hilbert space, one has the following theorem known as 
the Marchenko-Pastur theorem 0 . 

Theorem: Let us assume for simplicity that the compo¬ 
nents of Xi are Gaussian random variables, zero Mean, 
unit variance, and bounded moments that is there is some 
bound B, independent of m, such that\/n, Edxijl’^) < B. 
Then n depends on m in such a way that m/n ^ r < I 
as n —>■ oo. Under these assumptions, the distribution of 
the eigenvalues of Eyy^ asymptotically approaches the 
Marcenko-Pastur law os n —>■ oo 


fix) 


y/{x - a){b - x) 
2TTxr 


( 5 ) 


where a = (I — i/r)^ and 5 = (1 -I- v^)^. For the case 
r = 1 it reduces to a famous quarter circle law. 
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One important question that arises is how does one 
compare the covariance matrix constructed from the ran¬ 
dom matrices to a covarience matrix obtained from a 
definite state? First thing, if one looks at this two co- 
varince matrices they are Gaussian. Then the ques¬ 
tion of randomness in the quantum system, this comes 
from the wave function. This is shown by one of the 
authors Q, the connection between the action equa¬ 
tion for the stationary points of the random matrix en¬ 
semble V{x) — — xi\) where V{x) is the 

potential and the quantum momentum function [^, Q 
p = x-xk moving poles are sim¬ 

ple poles with residue —ih (we take here h = m = 1) 
and Q{x) is the singular part of the quantum momentum 
function. 

Therefore, if one takes this point of view that there is 
randomness associated with the wave function then one 
can see that the Wigner function or distribution is like 
the Wishart distribution. In random matrix theory the 
ensembles constructed using covariance matrix are the 
known as the Wishart Matrices and their joint proba¬ 
bility distribution function of eigenvalues is given by the 
Laguerre polynomials The Wigner function for the 
Harmonic oscillator are the Laguerre polynomials 0- 
It is shown in the reference m the Wigner function of 
the squeezed displaced vacuum state gives the Laguerre 
polynomials. 

It is well known now that the squeezed states are non 
classical states and shows the property of entanglement 
Q . The condition for squeezing is characterized in terms 
of covariance matrix usi ng Wig ner function reduces to a 
uncertainty relationship |12 h14| | 

v+'-n>o, (6) 

where V is the covariance matrix and is obtained by 
arranging the phase space variables and the hermitian 
canonical operators into four-dimensional column vectors 

C = ( 91 Pi 92 P2 ) , I = ( 9i Pi 92 P2 ) ■ 

Then the vector ^ satisfies the following commutation 
relationship (T^ 

[Caj ~ i Oi, P = 1, 2, 3,4; 

= (o 5)'•'=(- 1 S)- 

In particular, the condition for the multimode squeezing 
is also defined in terms of covariance matrix in ref [T 3 | . 

Let us consider a ensemble n-bipartite system of 
two modes described by annihilation operators = 

I]r=i(9^i + ipAi)/V2 and hi = I]"=i(9Bi + ipSi)/^- 
By definition, a quantum state p of the bipartite system 
is separable if and only if p can be expressed in the form 

n 

P= X! 

Ai,Bi = l j 


with nonnegative Pj’s, where PjAi’s and PjB^s are den¬ 
sity operators of the modes of A and B respectively. To 
check that the given n-bipartite system of two modes is 
entangled or not we would perform that partial transpose 
operation on the density matrix ([8]) denoted PT. Then 
one calculates the eigenvalues partial transposed density 
matrix. If all the eigenvalues are positive then the it is a 
separable density matrix if it has atleast one eigenvalue 
then the state is entangled. This is the Peres-Horodecki 
separability criterion. 

In the case of large matrices it is hard to check the pos¬ 
itivity of the all eigenvalues, but the Marchenko-Pastur 
theorem in random matrix theory gives a condition for 
the convergence of the eigenvalues for the large matrices. 

This property of the random matrices is used to check 
weather a given state is separable or entangled. 

In order to apply the Peres-Horodecki separability cri¬ 
terion for the Gaussian states one has to go over to phase 
space picture and then study the partial transpose oper¬ 
ation in the Wigner picture, it is convenient to arrange 
the phase space variables as 

C = ( 9A1 • ■ ■ 9Art 9B1 ■ ■ • gBn PAl ' ' ' PAn PBl ' ' ' PBn ) ■ 


C = ( QAl ■ ■ ■ qAn QBI ■ ■ ■ gBn PAl ' ' ' PAn PBl ' ' ' PBn ) ■ 
The commutation relations take the compact form[T^ 

[Cq; ^/3] “ i ^q/3 : /5 = 1, 2,3,4; 


/ J.. 

• 


■ • 1 T 0 1 \ 


■ ■ ’ (-1 oj 

Vo •• 

• J J 


here is a 4n x 4n matrix. In the two mode case it 
has been shown in ref under partial transpose opera¬ 
tion on the bipartite density operator transcribes has the 
following transformation on the Wigner distribution: 

PT : W{qi,pi,q2,P2) — W{qi,pi,q2,-p2)- (10) 

It is clear from the above that the sign momentum of 
the second is changing. Therefore when we apply partial 
transposition on multi-mode system is equivalent to flip¬ 
ping the sign of momentum variable of second mode that 
is PBi ■ 

The Peres-Horodecki separability criterion for Gaus¬ 
sian state is defined in ref Q as reads: if p is separable, 
then its Wigner distribution necessarily goes over into a 
Wigner distribution under the phase space mirror reflec¬ 
tion A. IT(A^), like IF(^), should possess the “Wigner 
quality”, for any separable bipartite state. 

In multimode case, the covariance matrix will be of 
dimension 4n x 4n that is dimension 2n from mode A 
and another dimension 2n from mode B. By following 
the procedure suggested in the ref Q and arranging the 
uncertainties or variances into a 4n x 4n real variance 


In terms of hermitian canonical operators 


( 8 ) 
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matrix V, defined through Vap = ({A^q,, A^^}), where 
A^ = ^ — ( £) , h ere (^q.) = tr^a/5. Then the uncertainty 
principle [3|, ll3L Il4| . 


V +-Vl>Q. 


( 11 ) 


Note that (6t implies, in particular, that V > Q, for de¬ 
tails ref d, 3, [iJI ■ 

Now by using a real linear transformation on ^ specified 
by 4n x 4n real matrix 


where A = (1, • • •, 1,1, —1, • • •, 1, —1), there are total 4 ?t, 
ones, the first 2n the -1-1 corresponds to for the qAi,QBi 
in the second 2n the -1-1 corresponds to the pAi and the 
— 1 corresponds to the pBi- In the covariance matrix this 
is equivalent to transposing each aij. Hence one has 


/-il ■■ 




\ ^nl ’ ’ 

^nn j 


(16) 


I ^ i' = (12) 


This transformation preserves the commutation relations 
M- We use in such a way that the 


C = { 9A1 <?B1 ■ ■ • qAn qsn PAI PBI ' ' ' PAn PBn ) • 

In terms of annihilation and creation operators which is 
equivalent to 

C' = ( OAl bBl ■■■ aAn bun VbI • • • ’ 

Thus the covariance matrix is written as 


V = 


/ CTll 


<^nl 



(13) 


where tTy is a 2 x 2 matrix constructed from expec¬ 
tation values of (aAi, bBi, aA^, b^^,) elements. Since V 
is a bonifide the covariance matrix of Gaussian states 
then the eigenvalue distribution of this covariance ma¬ 
trix should follow Merchnko Pasture law. By identifying 
this covariance matrix Gaussian states with the Gaus¬ 
sian random covariance matrix we get that m = 2n and 
n = 4n. One can clearly see that as m —>■ oo and n —> oo 
the ration m/n = 1/2. Therefore, the distribution of the 
eigenvalues oi V = asymptotically approaches the 

Marcenko-Pastur law as m, n —>■ oo 


f{x) 


y/{x - a){b - x) 

7TX 


(14) 


where a = (1 — y and 6 = (1 -h y ^)^. It is evident 
from eq (nni) that under the partial transpose the Wigner 
distribution undergoes mirror reflection, and that reflects 
at the covariance matrix level by the following transfor¬ 
mation V ^ V = APA. Then the new distribution func¬ 
tion 1T(A,J) has to be a Wigner distribution if the state 
under consideration is separable, then one has 


V+^i^>0, V = AVA, 


(15) 


Given a Gaussian states say p constructed form Hilbert 
space Ha ® Hb, under partial transposition of H, p 
has to be a bonifide density matrix. Then the Peres- 
Horodecki criterion is imposed on the covariance matrix 
constructed from the phase space variables Q . Then the 
Peres-Horodecki criterion for the multimode continuous 
variables or the Simon’s criterion for multimode reduces 
to the following theorem. 

Theorem : It is necessary and sufficient condition for 
separability for a given multimode Gaussian state if the 
eigenvalue distribution of the partially transposed covari¬ 
ance matrix V satisfies Marcenko-Pastur law . 

Proof: Let us assume that the partial transposed co- 
variance matrix P is a bonifide covariance matrix of 
Gaussian states and the eigenvalue distribution of this 
covariance matrix satisfy’s the Merchnko Pasture law. 
Thus by identifying this covariance matrix of Gaussian 
states with the Gaussian random covariance matrix we 
get that m = 2n and n = 4n. One can clearly see that as 
TO —>■ oo and n —>■ oo the ratio m/n = 1/2. Therefore, the 
distribution of the eigenvalues of P = asymptoti¬ 
cally approaches the Marcenko-Pastur law as to, n —>■ c» 

m = (17) 


where a = (1 — y 1)^ and 6 = (1 -f y 1)^. This gives a 
condition that the eigenvalues x is bounded above and 
below by 2^/2 < x — 3 < —2'j2. Within this range the 
states are separable and outside they are entangled. 

In Conclusion, we extend Simon’s criterion for Gaus¬ 
sian states to the multi-mode Gaussian states using the 
Marchenko-Pastur theorem from random matrix theory 
as necessary and sufficient condition for separability. 
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